function ejer4

	n = 3;

	b1 = [ones(n,1);zeros(n,1)];
	
	b2 = ones(2*n,1);
	
	I = eye(n);
	
	B = diag(4*ones(n,1))-diag(ones(n-1,1),1)-diag(ones(n-1,1),-1);
	
	fprintf('matriz A\n:')
	A=[B,-I;-I,B]
	
	fprintf('vector b1 y b2 \n:')
	b1 = [ones(n,1);zeros(n,1)]		
	b2 = ones(2*n,1)
	
	
	fprintf('solucion con sollu para b1: x\n')
	x1 = sollu(A,b1);
	x1'
		
	
	fprintf('solucion con soleg para b1: x\n')
	x2 = soleg(A,b1)
		
	
	fprintf('solucion con sollu para b2: x\n')
	x3 = sollu(A,b2);
	x3'
	
	
	fprintf('solucion con soleg para b2: x\n')
	x4 = soleg(A,b2)
	
	
endfunction
